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We demonstrate how the transverse t + t~ spin correlations can be used 
to determine whether a decaying Higgs boson is a mixed CV eigenstate, 
thereby directly probe the presence of CV violation in the neutral Higgs 
boson sector. We investigate the subsequent decay chain H — > t + t~ — > 
p + v T p~v T — > 7r + 7r°P r 7r _ 7rV r . The prospects for the measurement of the 
pseudoscalar admixture in the Ht + t~ coupling to a Standard Model Higgs 
boson with a mass of 120 GeV are quantified for the case of e + e~ collisions 
at 350 GeV center-of-mass energy and 1 ab~ l integrated luminosity. The 
Standard Model Higgsstrahlung production process is used as an example. 

PACS numbers: 14.60.Fg, 14.8d.Bn, 14.80. Cp 

1. Introduction 

A most distinctive feature of the extended models of Standard Model 
(SM) such as Minimal Supersymmetric Standard Model (AiSSM.), or gen- 
eral Two Higgs Boson Doublet Model {2TLVM.), is the existence of addi- 
tional Higgs bosons, i.e. a charged Higgs boson pair, , one CV odd scalar 
A and two CV even scalars h°,H°. Their mass and coupling patterns vary 
with the model parameters. The discovery of a neutral Higgs boson with 
mass in the range 114 GeV < Mh < 140 GeV will raise the question 
whether the observed particle is the SM Higgs boson or the lightest boson 
from the Higgs boson sector of a SM extension. Even if the one doublet 
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SAi turns out to be a good model, we will want to strictly determine that 
the observed Higgs boson is indeed CV even in its nature. 

Whichever scenario is realized in Nature, the major goal of future high 
energy experiments will be to distinguish among these models. The precision 
measurements of all relevant Higgs boson couplings are powerful to obtain 
information about additional Higgs doublets, their structure and masses. 
Such measurements can only be performed at the e + e _ linear collider (LO) 
such as TESLA. The determination of the quantum numbers of the Higgs 
boson(s) and sensitivity to the CV violation represent a crucial part of this 
test and will allow to establish the Higgs boson mechanism as the mechanism 
of the electroweak symmetry breaking. 

The determination of the parity and the parity mixing of spinless Higgs 
bosons have been extensively investigated in Refs. [1-8]. The model in- 
dependent identification of parity of the Higgs particle has recently been 
demonstrated for H/A° — > r + r + decay chain in Refs. [9-13]. 

In Ref. [9], see also Ref. [12] the reaction chain e + e~ — > Z(H/A°), 
H/A° — > t + t - , — > -K^Ur^Ur) was studied. It was found that even small 
effects of smearing seriously deteriorate the measurement resolution. How- 
ever, the spin effects of the decay chain H/A° — > t + t~ — > p + v T p~v T — ► 
7r + ir v T ir~ir°v T give a parity test independent of both model (e.g. SM, 
MSSM) and Higgs boson production mechanism (e.g. Higgsstrahlung, 
WW fusion). Using reasonable assumptions about the SM. production cross 
section and about the measurement resolutions we have found that, with 
500 fb^ 1 of luminosity at a 500 GeV e + e~ linear collider, the CV of a 120 
GeV Higgs boson can be measured to a confidence level greater than 95%, 
see Ref. [10]. 

In Ref. [11] we demonstrated that a measurement of the r impact param- 
eter in one-prong r decay is useful for the determination of the Higgs boson 
parity in the same decay chain. For a detection set-up such as TESLA, use 
of the information from the r impact parameter can improve the significance 
of the measurement of the parity of a Standard Model 120 GeV Higgs boson 
to ~ 4.5(7. 

In this paper we continue to investigate the case of a Higgs boson decay 
into r + r~ and we concentrate on the more general case where mixed scalar 
and pseudoscalar couplings of the Higgs boson to r leptons are simultane- 
ously allowed, see Ref. [13] for details study. 

The rest of the paper is organized as follows. In Section 2 we present 
general information about measurement of the CV quantum numbers of 
Higgs boson using rr decay channel. In section 3 we recall basic properties 
of the density matrix for the pair of r leptons produced in Higgs boson 
decay. In Section 4 we define an observable. Numerical results are presented 
in Section 5. The Summary closes the paper. 
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2. CV properties of the Higgs boson from tt decay 



Various more or less complicated extensions of the SM have been pro- 
posed in literature, see e.g. Ref. [14-20], but their common feature is that in 
case of models with CV violation in the Higgs boson sector the mass eigen- 
states of the neutral Higgs bosons are not precisely CV eigenstates h°, H° 
and A . CV violation results in three neutral states of mixed CV character. 
For these models to estimate precision of the H — A mixing angle is essen- 
tial. We consider a potential of establishing CV properties of Higgs bosons, 
without any assumption about a CV violation model, from the analysis of 
the angular distributions of the r + r~ decay products in the plane trans- 
verse to the t + t~ axes. The transverse spin effects in r pair production, 
reflected in correlations between r decay products, are helpful to distin- 
guish between the scalar J vc = ++ , pseudoscalar J vc = h and mixing 
natures of the spin zero (Higgs) particles. When the Higgs boson is light 
enough that the W + W~ decay channel remains closed, the most promising 
decay channel of SM neutral Higgs particle sensitive for spin correlations, 
is the t + t~ mode. The t + t~ channel is useful in the SM for Higgs masses 
less than ~ 140 GeV. Up to this mass, the Higgs particle is very narrow, 
T(Hsm) < 10 MeV with BTZ(H S m r+r") ~ 9%. In Supersymmetric 
theories, the t + t~ channel is useful over a much larger mass range. In this 
approach the CV properties of Higgs boson can be studied independently of 
a production mechanism and the specific model and can be considered as 
the most general one. 

It is generally believed that a Monte Carlo simulation of the full chain 
from the beam collision to detector response is the most convenient tech- 
nique to investigate such studies. In our analysis, we will take as an example 
the e + e~ — ► ZH; Z — ► p + p~; H — > r+r~ production process. We discuss 
a method for the CV quantum numbers measurement of the Higgs boson 
with a mass of 120 GeV for the case of e + e~ collisions with yfs = 350 GeV 
using H/A° — > t + t~; — > p ± i> T (u T ); /r* 1 — ► tt^tt decay chain. All the 
Monte Carlo samples have been generated with the TAU0LA Monte Carlo 
library [21-23]. For the production of the r lepton pairs the Monte Carlo 
program PYTHIA 6 . 1 is used [24]. The effects of initial state bremsstrahlung 
were included in the PYTHIA generation. For the r lepton pair decay with 
full spin effects included in the H — ► r + r~; — ► p^v T {y T )\ — > Tr^-zr 
chain, the interface explained in Refs. [9, 10] was used. It is an extended 
version of the standard universal interface presented in Ref. [25], see also 
Ref. [26]. 
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3. Mixed scalar pseudoscalar case in Monte Carlo algorithm 

Let us now, only very briefly repeat the basic information about the spin 
correlations and their implementation in our Monte Carlo algorithm. The 
detailed description of the method can be found in Ref. [21] and the full 
description of the algorithm is given in Ref. [9]. 

The main spin weight of our algorithm for generating the physical pro- 
cess of r lepton pair production in Higgs boson decay, with subsequent decay 
of r leptons as well, is given by 

wt 4 1 + EEM^]. (J) 

V i=ij=i 

where hi and hi are the polarimeter vectors that depend respectively on 
r ± decay products momenta; R^ is the spin density matrix. For the mixed 
scalar-pseudoscalar case, when the general Higgs boson Yukawa coupling to 
the r lepton 

f(a + i&7 5 )r (2) 
is assumed, we get the following non-zero components of Rij, see Ref. [13]: 

a 2 [3 2 -b 2 2abp . , 

i? 33 - -1, flu - R22 - a2(j2 - b2 , R12 - -R21 - q2(32 - b2 , (3) 



where j3 = \Jl — Am 2 /m 2 H . The crucial point is that, in general, a and b 
are of comparable magnitude in a CP violating extension of S M . For a CP 
conserving Higgs sector, either a = or b = implying a pure pseudoscalar 
or scalar case respectively. For a CP mixed eigenstate, both a and b are non 
zero. Thus any significant deviation of R\2 or equivalently R21 from zero 
provides a signature for CP violation in the Higgs sector independent of the 
specific model. If we express Eq. (2) with the help of the scalar-pseudoscalar 
mixing angle (f>: 

f N (cos (f> + i sin (^75 ) r, (4) 

the components of the spin density matrix can be expressed in the following 
way: 

cos 4> 2 1 - sin 4> 2 2cos0sin0/3 

^11 — — To a2 , : To" 1 K 12 — --K21 — Jo Q2 1 : — If" ^1 

cos (p z p z + sin cp z cos (fi z p z + sin <p z 

In the limit (3 — > 1 these expressions reduce to the components of the rota- 
tion matrix for the rotation around the z axis by an angle — 2(j>: 



R11 = R22 = cos 2(f>, R\ 2 = —R21 = sin 20. 



(6) 
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4. Definition of the observable 

Let us now recall a observable which we have introduced to distinguish 
between scalar-pseudoscalar mixed state of the Higgs boson. The method 
relies on measuring the acoplanarity angle of the two planes, spanned on 
decay products and defined in the p + p~ pair rest frame. For that purpose 
the four-momenta of and 7r° need to be reconstructed and, combined, 
they will yield the p^ 1 four-momenta. All reconstructed four-momenta are 
then boosted into the p + p~ pair rest frame. The acoplanarity angle p*, 
between the planes of the p + and p~ decay products is defined in this 
frame. The correlation, in the case of the Higgs boson of combined scalar 
and pseudoscalar couplings with the mixing angle <f>, is between transverse 
components of t + spin polarization vector and transverse components of 
r~ polarization vector rotated by an angle 2(j). Therefore the full range of 
the variable < p* < 2tt is of physical interest. To distinguish between the 
two cases <p* and 2tt — p* it is sufficient, for example, to find the sign of 
p n - • n + , where n + is a vector normal to the plane spanned by the visible 
decay products of p + , n + = p n + x p^o . The range < p* < it corresponds 
to the negative sign case, otherwise one should make the replacement p* — > 
2ir — p* . If no separation was made, the parity effect, in case of mixed Htt 
coupling, would wash itself out. Additional selection cuts need to be applied. 
Otherwise the acoplanarity distribution is not sensitive to transverse spin 
effects at all. The events need to be divided into two classes, depending on 
the sign of yit/2, where 

E n + + E^o E w - + E^o 

The energies of vr^,^ are to be taken in the respective r ± rest frames. 
In Refs. [10, 11] the methods of reconstruction of the replacement rest 
frames were proposed with and without the help of the r impact parameter. 
We will use these methods here as well, without any modification. 



5. Numerical results 

Let us turn our attention to the measurable distributions. We have used 
the scalar-pseudoscalar mixing angle <j) = f and, as the reference, we have 
used the pure scalar case = 0. In our study, that is for the 350 GeV e + e~ 
center-of-mass energy, Higgs boson mass of 120 GeV and Higgsstrahlung 
production process, we took N a = 62.7 • 10~ 3 [fb] for the scale of the plot. 
However, in general case 

N a = ^-atotai{e + e- XH)BTZ{H t + t-){BTZ{t -► pu T )f (8) 
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Fig. 1. The acoplanarity distribution (angle ip* ) of the p + p~ decay products in the 
rest frame of the p + p~ pair. Gaussian smearing of n 's momenta and generator 
level rest frames are used. The thick line corresponds to a scalar Higgs boson, 
the thin line to a mixed one. The left figure contains events with y.iy.2 > 0, the 
right one is for 1/1J/2 < 0. 



is a suitable choice. 

In Fig. 1 the acoplanarity distribution angle ip* of the p + p~ decay prod- 
ucts which was defined in the rest frame of the reconstructed p + p~ pair is 
shown. A detector-like set-up is included in exactly the same proportion as 
in Ref. [13]. Unobservable generator-level rest frames are used for the 
calculation of selection cuts. The two plots represent events selected by the 
differences of tt^t: energies, defined in their respective rest frames. In 
the left plot, it is required that yiU2 > 0, whereas in the right one, events 
with yiy2 < are taken. This figure quantifies the size of the parity effect 
in an idealized condition, which we will attempt to approach with realistic 
ones. 

The size of the effect was substantially diminished when a detector-like 
set-up was included for r ± rest frames reconstruction as well, see Fig. 2. 
The general shape of the distributions however remained. At the cost of 
introducing cuts, and thus reducing the number of accepted events, we could 
achieve some improvement of the method, as in Ref. [11]. If we require the 
signs of the reconstructed energy differences y\ and 1/2, see Eq. (7), to be the 
same whether the method is used with or without the help of the r lepton 
impact parameter, only ~ 52% of events are accepted. The relative size of 
the parity effect increases. Results are presented in Fig. 2. 

Both distributions clearly distinguish between decays of scalar or mixed 
Higgs boson. From the measurement of these distributions we can establish 
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Fig. 2. The acoplanarity distribution (angle tp* ) of the p + p~ decay products in the 
rest frame of the p + p~ pair. Gaussian smearing of it's and Higgs boson momenta, 
are included. Only events where the signs of the energy differences y\ and yi are 
the same, if calculated using the method without or with the help of the r impact 
parameter are taken. The thick line corresponds to a scalar Higgs boson, the thin 
line to a mixed one. The left figure contains events with y\yi > 0, the right one is 
for yiy 2 < 0. 

the CV properties of the Higgs boson. More precisely, we have found, see [13] 
for details, that for an integrated luminosity of 1 ab~ 1 , at 350 GeV center- 
of-mass energy, a high precision LC detector such as the proposed TESLA, 
should be able to measure the scalar-pseudoscalar mixing angle for the Htt 
coupling with 6° accuracy in the case of a Standard Model Higgs boson mass 
of 120 GeV. 

Summary 

We have studied measurement opportunities of the CV properties of a 
120 GeV Higgs boson at an e + e~ collider, e.g. TESLA. Our results show 
that if the Higgs sector is CV violating then there is a possibility to explicitly 
test this CV violation trough spin correlations between final state particles 
in H — ► t + t~; r — > p v T (v T ); — > tt^tt decay chains. We have found 
that the mixing scalar-pseudoscalar angle can be determined with statistical 
precision of 6° for the case of e + e~ collisions of 350 GeV center-of-mass 
energy with an integrated luminosity of 1 ab~ l and for Higgs boson mass 
of 120 GeV. However, we assume that the e + e~ — > ZH; Z — > p + p~; 
H — > t + t~; — > p~^v T (y T ); p^ — > vr^vr decay chain is background free. 
We have not introduced any cuts etc. that might be required to guarantee 
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this. 

Finally, let us note that this method can be applied to measure the parity 
properties of other scalar particles, not necessarily only Higgs boson(s). 
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